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ABSTRACT 

Discerning the exact nature of the sub-mJy radio population has been historically diffi- 
cult due to the low luminosity of these sources at most wavelengths. Using deep ground 
based optical follow-up and observations from the Spitzer Space Telescope we are able 
to disentangle the radio-selected Active Galactic Nuclei (AGN) and Star Forming 
Galaxy (SFG) populations for the first time in a deep multi-frequency VL A/MERLIN 
Survey of the 13 H XMM-Newton/ Chandra Deep Field. The discrimination diagnostics 
include radio morphology, radio spectral index, radio/near-IR and mid-IR/radio flux 
density ratios. We are now able to calculate the extragalactic Euclidean normalised 
source counts separately for AGN and SFGs. We find that while SFGs dominate at the 
faintest flux densities and account for the majority of the up-turn in the counts, AGN 
still make up around one quarter of the counts at ~ 50 /zJy (1.4 GHz). Using radio 
luminosity as an unobscured star formation rate (SFR) measure we are then able to 
examine the comoving SFR density of the Universe up to z = 3 which agrees well with 
measures at other wavelengths. We find a rough correlation of SFR with stellar mass 
for both the sample presented here and a sample of local radio-selected SFGs from the 
6df-NVSS survey. This work also confirms the existence of, and provides alternative 
evidence for, the evolution of distribution of star formation by galaxy mass: "down- 
sizing". As both these samples are SFR-selected, this result suggests that there is a 
maximum SFR for a given galaxy that depends linearly on its stellar mass. The low 
"characteristic times" (inverse specific SFR) of the SFGs in our sample arc similar to 
those of the 6dF-NVSS sample, implying that most of these sources are in a current 
phase of enhanced star formation. 
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1 INTRODUCTION 

Whilst observations of powerful, radio-loud Active Galac- 
tic Nuclei (AGN) were an early probe of the distant Uni- 
verse fsee lStern fc Spinradlll999l . for a summary), starburst 
galaxies are three or more orders of magnitudes less lumi- 
nous at radio wavelengths and hence difficult to observe 
at large distances. However the deepest r adio surveys at 
1.4 GHz now reach an rms below 10 /iJy l|Richardsl I2OO0I: 
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Seymour et ail |2004 iBiggs fc Ivisonl 120061 : iFomalont et ail 
200rj . etc.). These deep radio observations reveal a (very 



well characterised) up-turn in the Euclidean normalised 
sources counts below 1 mjy above that predicted from the 
extrapolation of the AGN counts measured at brighter 
flux densities. This up-turn has been attributed to the 
emergence of a star forming galaxy (SFG) population, 
requiring strong evolution of the SFG radio luminosity 
function ([R owan-Robinson e t al. 



ISevmour et al J [20041 ; iMoss et~ai 



le Q- r vj radio l uminosity 
19931: IHopkins et alJll998l: 



20071 ) . although some au- 



thors argue that there is a significant contribution due to rel- 
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atively weak radio AGN l|Simpson et al.ll2006l ; llluvnh et al.l 
l2007l : lBarger et al.ll2007h . 

Determining the nature of individual radio sources has 
remained difficult due to their low luminosities at other 
wavelengths. Discrimination between radio-selected SFGs 
and AGN in these deep surveys would, for example, al- 
low an independent measure of the star formation history 
of the Universe as star formation rate (SFR) is directly re- 
lated to radio l uminosity for g alaxies with weak or no AGN 
radio emission (Condon 1992). The comoving SFR density 
(SFRD) has been de termined previously from deep radio 
data out to z ~ 1.6 (|Haarsma et al.l 120001 ). but the study 
presented here offers several improvements and advantages, 
noticeably a more robust and systematic discrimination be- 
tween AGN and SFG, the use of more accurate redshifts, an 
extension to higher redshift and a larger area/sample size 
by a factor of ~ 5. 

Whilst radio surveys have not traditionally been used to 
determine SFRs in the distant Universe, other methods have 
their limitations. The far-IR and sub-mm wavelength range 
is limited by the relatively poor sensitivity and angular res- 
olution of current instrumentation and telescopes. The mid- 
IR, whilst very sensitive, is not as reliable a tracer of SFR at 
high redshift as the far-IR and sub-mm (e.g. lPapovich et al.l 
120071 ). UV emission from young massive stars can easily be 
obscured by dust, but the extrapolation of models derived 
from local galaxies likely break down for the more lumi- 
nous galaxies, > L», found in deep optical surveys. Fur- 
thermore, the rest-frame UV is redshifted across a large 
wavelength range (UV/optical/near-IR) at cosmological dis- 
tances so different detectors are needed to trace the UV at 
different redshifts. Determination of the SFR from optical 
emission lines is also affected by obscuration and shifting 
to longer wavelengths and hence in and out of observable 
wavelength windows. The optical emission lines can also suf- 
fer from absorption by dust. X-ray emission associated with 
star formation (e.g. from binary stars) is intrinsically weak 
and is subject to varying degrees of photoelectric absorp- 
tion. Even at the faintest X-ray fluxes reached in the CDFN 
(~ 10 ~ 17 ergcm -2 ^" 1 ), the source counts are dominated by 
AGN l|Bauer et al.| [2004). Given that forthcoming radio fa- 
cilities (e.g EVLA, e-MERLIN, LOFAR, ASKAP, and even- 
tually the SKA) will reach many orders of magnitude deeper 
than the current deepest radio survey, and will be able to 
detect star forming galaxies out to the era of re-ionization 
(z < 7), it is timely to consider how to characterise the emis- 
sion from faint radio sources and which supporting data are 
the most valuable. 

We make an important distinction here between radio 
indicators of an AGN and non-radio indicators of an AGN. 
We need to separate the radio sources into those whose radio 
emission is AGN dominated and those that are consistent 
with being dominated by star formation. By AGN powered 
indicators we mean some measure that indicates that the ra- 
dio emission from a galaxy is dominated by accretion onto 
a super-massive black hole and/or from associated jets and 
lobes (i.e. radio-loud AGN). We assume, for simplicity, that 
the dominant power in most sources is either from an AGN 
or from star formation, but otherwise make no a priori as- 
sumptions about the sources, e.g. the SED. The radio AGN 
indicators we use in this paper are radio/near-IR flux den- 
sity ratio, mid-IR/radio flux density ratio, radio morphol- 



ogy and radio spectral index. We do not include non-radio 
AGN indicators, e.g. classical methods like optical/IR emis- 
sion lines, optical morphology, X-ray observations etc. and 
methods from recent results with the Spitzer Space Telescope 
based on mid-IR SEDs (e.g. iDonlev et aill2005l ; [stern et all 
2005| lLacv et alj|2003 ) . As these non-radio indicators tell us 
nothing directly about the nature of the radio emission we 
do not use them in our discrimination meth ods. A compa- 
rable philosophy was previously adopted by iMuxlow et all 
|2005l ). 

In this paper we will present our method for classifying 
the faint radio sources and examine the properties of the 
SFG population. The radio, optical and IR observations are 
summarised in section 2. The various methods of discrimi- 
nating between AGN and SFGs are presented in section 3. In 
section 4 we present the Euclidean normalised source counts 
separately for AGN and SFGs for the first time. In section 
5 we study the general properties of the radio-selected SFG 
population, present the comoving SFR density of the Uni- 
verse up to redshift ~ 3 derived from our data, examine 
the distribution of star formation with host galaxy stellar 
mass and derive characteristic times. We conclude this ar- 
ticle in section 6. Throughout we use a concordance model 
of Universe expansion, 11 m = 1 — Oa = 0.3, Qo = 1, and 
H = 70 kms^Mpc -1 (jSpergel et al.ll2003h . Magnitudes 
are AB unless otherwise stated. 



2 THE DATA 

The 13 H XMM-Newton/ 'Chandra Deep Survey field presents 
us with a unique data set for performing a radio-based 
SFG/ AGN separation of radio sources. This field, centred 
at 13 h 34 m 37 s +37° 54' 44 " was the location of one of the 
deepest ROSAT surveys (|M C Hardy et all Il998l ). and lies 
in a region of extremely low Galactic absorption (Nh ~ 
1 cm" 



2.1 The radio data 

2.1.1 The VLA 14 GHz observations 

The 13 H XMM-Newton/ Chandra deep field was observed by 
the VLA at 1.4 GHz for 14hrs in A-array in August 1996 
and 10 hrs in B-array in November 1995. A catalogue of 449 
sources above 30 fiJy (4 a detection limit) in a 30 arcmin di- 
ameter circular field of view (0.196 deg 2 ) was obtained with 
a resolution of 3.3 arcsec. The Euclidean normalised source 
counts, corrected for the incompleteness o f the catalogue 
due to instrumental effects, were presented in lSevmour et all 
(|2004 hereafter S04). The 449 sources detected at 1.4 GHz 
with the VLA constitute the parent sample on which this 
work is based. In contrast to S04 we use the peak flux density 
instead of total flux density for sources that are significantly 
resolved and at low S/N as the peak flux is a better measure 
of the true flux in these cases. 

The size of this sample (in terms of area and number of 
sources) is well positioned between those fields of the sim- 
ilar area but are slightly deeper (by ~ 40% e .g. the Lock- 
man Hole and the HDFN. lBiggs fc Ivisonll2006l) and those of 
larger area (e.g AEGI S and COSMOS. ITvison et al.ll2007al ; 
ISchinnerer et al.ll2007l ), but which are less deep by ~ 40%. 
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2.1.2 The MERLIN 1.4 GHz observations 

The 13 H field was observed with Mult i- Element Radio- 
Linked Network in April 1999. Due to the smaller field of 
view, 10 arcmin, four MERLIN pointings were used to cover 
most of the VLA field of view. MERLIN has a higher resolu- 
tion, 0.2 arcsec, than the VLA due to its greater maximum 
base-line. We used the parent catalogue as a detection list 
and only imaged the positions around known sources (in- 
cluding bright sources not in the principle VLA or MER- 
LIN field of view). We made MERLIN images for all VLA 
detections with Si.4GHz > 72 /xjy that lay inside the MER- 
LIN coverage. These images were combined by first mak- 
ing maps without the beam deconvolved and beams around 
each source in the sky plane, then averaging the maps and 
beams before deconvolving the average beam from the av- 
erage map as described by A. Zoghbi et al. (in prep.) The 
combination of datasets from different telescope arrays is 
non-trivial as the four MERLIN pointings are offset from 
the VLA phase center and both the VLA and MERLIN im- 
ages suffer from radial smearing. However for sources in the 
parent sample having flux densities above ~ 100 /iJy, and 
that are favourably positioned, we can make high-resolution 
maps from combined MERLIN and VLA data. 



2.1.3 The VLA 4 . 8 GHz observations 

The 13 H field was observed with the VLA at 4.8 GHz in 
April 1991. The observations cover the whole field of view 
of the 1.4 GHz data, 30 arcmin (observatio ns originally de- 
signe d to match the ROSAT field of view lM c Hardv et all 
1 19981 ) with 51 separate pointings. The observations reach 
33 /iJy rms, and were performed in VLA D-array hence the 
restoring beam is 14 arcsec, around a factor four gre ater than 
the 1. 4 GHz observations. From the catalogue of ISevmourl 
(2002) we find 45 sources with flux densities > 4 a and 
within the 30 arcmin diameter field of view of the 1.4 GHz 
observations. Of these sources 45 have counterparts in the 
1.4 GHz catalogue within 7 arcsec and hence we can cal- 
culate their 1.4 — 4.8 GHz spectral indices. We define the 
radio spectral index, a, by S v oc v a , and hence a — 
— 1.85 log 10 (S , i.4/S , 4.86). We use total, rather than peak, flux 
densities when calculating the radio spectral index. In our 
earlier study (S04), we made tapered maps of the 1.4 GHz 
data but we did not find significant extended radio emission 
that could affect the calculated spectral indices. There is one 
instance where two nearby 1.4 GHz sources (separated by 
20 arcsec) are matched to a single 4.8 GHz source and there- 
fore we cannot calculate the spectral indexes separately for 
these two sources. 



2.2 Optical, near-IR and mid-IR data 

2.2.1 Optical and near-IR data 

We have obtained imaging in many bands from the 
near-UV to near-IR over recent years: u* , g' and i'- 
band from MegaCam/CFHT, B, R, I and z'-band from 
SuprimeCam/Subaru, Z-band from WFC/INT, J-band 
from WIRC/Palomar, #-band from WIRCam/CFHT, and 
JVband from WFCAM/UKIRT. This imaging covers the 
entire parent sample bar the WIRC J-band (which covers 



approximately 3/4 of the survey area), and has been sup- 
plemented by hundreds of optical spectra of the counter- 
parts to our radio and X-ray sources. Of our parent sample 
164/449 currently have redshifts determined from optical 
spectroscopy. 



2.2.2 Mid-IR data 

The 13 H field was observ ed by the IRAC dFazio et al]|2004f) 



and MIPS instruments (Ricke et al. 2004) on board the 



SpUzer Spaee Telescope Iwerner et al ' B in July 2005 



as part of MIPS instrument team GTO time (PI G. Rieke, 
program identification number 81). A strip of ~ 0.5 x 1 deg 
containing the 13 H field was imaged by both instruments. 
The IRAC data consisted of observations with all four chan- 
nels (3.6, 4.5, 5.6 and 8.0 /im) reaching 1 a flux density limits 
of 0.7, 1.2, 7.2 and 7.6 /iJy respectively. These values rep- 
resent the limiting flux density of a point source for which 
we can make a flux measurement in a 3.8 arcsec diameter 
aperture with > la accuracy. Most, 388/449 (87%), of our 
parent sample were detected in at least the 3.6 /im channel. 
The 24 /im data was taken as part of a simultaneous scan 
map in each of the three bands (24, 70 and 160 /im) reaching 
a 3<r detection limit of 98 /iJy in a 10.5 arcsec aperture. We 
found 330/449 (73%) of our parent sample had 24 /im coun- 
terparts. We do not use the 70 and 160 /im data here as they 
are not sensitive enough to provide information on most of 
the radio sources, and then only those at low redshift. 



2.2.3 Photometric Redshifts 

We have used our multi-band near-UV to near-IR together 
with IRAC 3.6, 4.5 and 5.8 /im data to calculate photometric 
redshifts for the optical counterparts to our radio sources. 
Full details of this process will be presented in Dwelly et 
al. (in prep, hereafter D08), but here we give a short sum- 
mary. Initial catalogues were generated independently for 
each optical, near-IR and m id-IR waveband using Sextrac- 
tor (|Bertin fc Arnoutsl 1 1996 1. These initial catalogues are 
then combined into a single "master" catalogue containing 
only unique detections. Aperture photometry is carried out 
in each waveband at the locations of each master catalogue 
source. In the optical/near-IR bands we used a 3 arcsec di- 
ameter aperture for bright sources (i' < 20), and a 2 arcsec 
aperture for fainter sources to improve S/N. In the IRAC 
bands a 5.8 arcsec diameter aperture was employed. Ap- 
propriate aperture corrections are applied gi ven the differ- 
ent p oint spread functions in each waveband (|Gawiser et alj 
120061 ). We match optical/near-IR counterparts to the parent 
radio sample are by searching for all objects in the "mas- 
ter" catalogue lying within 1.5 arcsec of the peak of the radio 
emission. Where available we use the high resolution images 
from combined MERLIN+VLA maps to determine the cor- 
rect optical counterparts. We examine the optical and ra- 
dio images for sources with more than one counterpart, and 
manually choose the most appropria te object. We then used 
the publicly available Hyperz code (|Bolzonella et all l200dh 
to determine photometric redshifts for each source from their 
multi-band aperture photometry. Only the wavebands short- 
ward of 5.8 fim were used to determine photometric redshifts 
in order to avoid complications with PAH features. We used 
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the standard set of synthetic evolving galaxy templates sup- 
plied with Hyper z. To ensure reliability, we ignore the pho- 
tometric redshifts calculated for all radio sources detected 
with > 3 a significance in fewer than 4 wavebands. Hence, 
60 radio sources from our parent sample have no photomet- 
ric redshift information including 31 with no counterparts 
detected in any optical, near-IR or mid-IR waveband. Most 
of the other 29 sources are detected only in the K— band 
data and/or the two shorter wavelength IRAC channels, and 
hence the derived photometric redshifts have very large un- 
certainties. We consider the likely nature of these objects in 
section 3.6.1. 

2.3 Local 6dF-NVSS comparison sample 

To illustrate our discrimination techniques we use the 6dF- 
NVSS sampl e of local spect r oscop ically identified radio 
sources from iMauch fc Sadler] (|2007T ) . These authors iden- 
tified ~ 8000 radio s ources from th e 1.4 GHz NRAO VLA 
Sky Survey (NVSS, ICondonl Il992l ) in the Second Incre- 
mental Data Release of the 6 degree Field Galaxy Sur- 
vey (6dFGSDR2) with galaxies brighter than K = 12.75 
(Vega) from the 2MASS Extended Source Catalogue. This 
sample covers about 17% of the sky and i nclude s galaxies 
across 0.003 < z < 0.3. IMauch fc Sadlerl (|2007f l discrimi- 
nated between AGN and SFGs using emission and absorp- 
tion features in the 6dF spectra and derive d local luminos- 
ity fun ctions for both populations using the Sa unders et al.l 
(1990) form of the luminosity function (LF) which is com- 
monly used at long- wavelengths (i.e. mid-IR and longward). 
Their sample has also been crossed-matched with the IRAS 
Faint Source Catalogue providing fluxes at 12, 25, 60 and 
100 /i m which they us e to co nfirm the radio-IR correlation. 
The IMauch k, Sadlerl (|2007l ) catalogue provides us with a 
local reference sample of SFGs with known radio, IR and 
K— band flux densities (from which we can derive SFRs 
and stellar masses) to which we can compare our results 
from higher redshift. 



3 RADIO-BASED EMISSION DIAGNOSTICS 

Here we present four radio emission diagnostics which di- 
rectly probe the physical origin of the radio luminosity: two 
of these are purely from the radio data (radio spectral in- 
dices, morphology), and two involve comparisons with ob- 
served flux densities at other wavelengths (mid-IR and near- 
IR) as a function of redshift. These diagnostics are mostly 
positive AGN discriminators, i.e. they generally imply that 
the radio emission of a source is very likely powered by AGN 
emission and not star formation. We do not include the ob- 
vious radio luminosity discriminator, where a radio lumi- 
nosity > 10 25 WHz" 1 would imply an unphysical SFR of 
~ 5000 M0yr -1 , as those eight sources clearly above this 
luminosity are flagged as AGN by one or more of the other 
methods below and potentially such objects do exist in small 
numbers at high redshift. 

We give the number of sources from our parent sample 
for which each of our AGN/SFG discriminators is applicable 
in Table [T] along with a short comment of why they apply 
to only such a number. We discuss the implications of these 
selection functions in t]3.5l 



Table 1. Number of sources from our parent sample for which 
each SFG/AGN discrimination method is applicable. We list dis- 
crimination method, the number of sources which it is applicable 
to and the reason why. 



Discrimination 


# 


Limitation of method 


morphology 


127 


brightest sources only and coverage 


spectral index 


45 


brightest sources only 


•524/^m/<S'l.4GHz 


445 


limited by 24 (im coverage 


S'l.4GHz/S , 2.2/im 


449 


applies to all sources 



3.1 Radio morphology 

Only six sources from the parent sample show the classical 
unambiguou s AGN signature of a double lobed, FR2 mor- 
phology (i.e. lFanaroff fc Rilevl[l974T ) in our 3.3 arcsec resolu- 
tion imaging. For other extended VLA sources there is not 
enough detail to distinguish between a clear AGN morphol- 
ogy and extended, galaxy scale emission from star formation. 

Additionally we include results of the combined 
MERLIN+VLA(A-array) maps at 1.4 GHz (A. Zoghbi et al., 
in prep.) which have a resolution of 0.5 arcsec. These results 
are limited to the brighter objects from our master sam- 
ple due to the sensitivity of our MERLIN images, and also 
the inner ~ 20 arcmins covered by the four MERLIN point- 
ings. We are therefore able to make combined maps for 127 
sources. Many show clear AGN jet/lobe morphology includ- 
ing the six sources flagged as AGN from the lower resolution 
VLA images. The other sources show either a very compact, 
unresolved morphology or are extended on ~ 1 arcsec scales, 
often almost completely resolved out due to the lack of sen- 
sitivity of MERLIN to extended structure. Although the un- 
resolved sources do not have significantly constrained bright- 
ness temperatures ( ^ 10 3 K) this compa ctness is most likely 
due to an AGN fe.g. lMuxlow et al]|2005l ). but a small chance 
remains that they could be due to nuclear starbursts. The 
sources with extended emission over the galaxy which is of- 
ten completely resolved out must be from star formation oc- 
curring throughout the galaxy. Hence, the MERLIN+VLA 
images show that 42/127 sources are very likely to be AGN. 

3.2 Radio spectral indices 

Using our 4.8 GHz data |Sevmourl f2002') we have radio spec- 
tral indices for 45 sources in our sample. These are among 
the brightest sources in our sample as the 4.8 GHz rms 
of 30 pJy corresponds to ~ 70 /xJy at 1.4 GHz assuming 
a*;! = —0.8. Radio emission from star formation, due to 
synchrotron radiation, will generally have a spectrum with 
a slop e in the range —0.5 < a\'% < —1-0 (|Thonrpson et all 
[2006J). Radio sources with spectral indices steeper than 
a = — 1 are typically targeted as high redshift radio galaxies. 
Therefore any radio sources with spectral indices not in the 
range -0.4 < atl < -1.1 (allowing for the errors in the flux 
densities) we flag as having AGN powered radio emission. 
We find 31 sources have spectral indexes not in this range 
(i.e. with atl > -0.4 or atl «S -1.1). 

This diagnostic is partic ularly powerfu l as it i s relatively 
insensitive to redshift; e.g. iKlamer et al.l l|2006l ) find that 
90% of a sample of high redshift AGN have radio spectra 



The Star Formation History of the Universe as Revealed by Deep Radio Observations 5 



best characterised as a power-law with no spectral steepen- 
ing. It is possible that high redshift SFGs could have flat- 
ter radio spectra due to free-free processes dominating the 
higher rest-frame frequencies, especially if the source is very 
young. However the 1.4 GHz counterparts to the 4.8 GHz 
sources are typically bright (> 1 mjy at 1.4 GHz) hence 
unlikely to include many high redshift SFGs. Most radio 
sources classified as AGN by their radio spectral index are 
also classified by other diagnostics (see Table |3j) . 



template tracks. Combining these uncertainties we get a 5 <r 
cut of 0.6 dex (= 5 x ^f(0.2 2 + 0.2 2 )/Zn(10)). 

We note that all bar one of the sources with unknown 
redshifts are well above our AGN/SFG discrimination line 
at almost all redshifts and hence are most likely all obscured 
AGN at high redshift (see § 13.5.11 for further discussion of 
these sources). 



3.3 Radio to near-IR flux density ratios 

"Radio-loudness" i s another classical dis criminator of ra- 
dio luminous AGN (|Kellerman et alii 19891 ) and traditionally 
based on radio to optical (observed or rest-frame) flux den- 
sity ratio. More recently radio-loudness has been defined in 
terms of absolute radio luminosity (e.g. iMiller et al.l 1 1990l . 
Li.4GHz > 10 25 WHz" 1 ), but this diagnostic would only ap- 
ply to eight of our parent sample and those have already 
been flagged as an AGN by one or both of the previous 
two diagnostics. The original definition of radio loudness 
was based on B-band observations, but given that at the 
redshifts in question, ^ z ^ 4, this band quickly starts 
sampling below the 4000 A break where variations in star 
formation rate and absorption by dust have a strong effect, 
we chose to use a longer wavelength band. In fact the only 
band short-ward of the mid-IR which remains longward of 
the 4000 A break up to z ~ 4 is K— band. Our A— band 
data are particularly deep: K — 22.6 (3a limit). 

In keeping with our philosophy stated in the introduc- 
tion, where we make no assumptions about the observed 
SED, we plot the observed flux density ratios against red- 
shift. Given that we do not know the nature of each radio 
source beforehand, we can only plot the un-k-corrected flux 
density ratios and compare them with tracks of different 
AGN and starburst galaxies. This approach avoids making 
any a priori assumptions about the nature of each source. In 
Figure[T](a) we present the 1.4 GHz to K— band flux density 
ratio against redshift. We plot the tracks of two star forming 
galaxies which have no strong AGN component. These are 
LIRG and ULIRG templates made from average composites 
from fits to about a dozen LIRGs and ULIRGs with very 
high fidelity observations. From 5 — 35 /im they are based on 
Spitzer IRS long low spectra. For 1 — 5 /im, the templates are 
built from stellar population models, with overall slopes con- 
strained by large beam IRAC and 2MASS ph otometry. The 
radio regime was determined with data from lCondon et"ahl 
|l99ll ). The corresponding SFRs of these two templates are 
approximately ~ 10 and ~ IOOMqJT - respectively. We 
also plot the tracks of classic al 'radio-loud' and 'radio-quiet' 
QSOs from the templates of lElvis et all (|l994h . 

Figure [1] (a) shows several clear radio-loud AGNs with 
flux density ratios similar to or greater than the radio-loud 
QSO template track. We can use this plot to define a lo- 
cus where a radio source will unambiguously identified as 
an AGN. We decide to take a cut 5 a above the higher SFG 
track (the ULIRG one) where a comes from taking an un- 
certainty of 20% in the radio flux density and 20% in the 
A'— band magnitude. We emphasise that these uncertainties 
are extremely conservative and that the majority of radio 
sources will have much smaller uncertainties, but we are 
also allowing for the uncertainty in the SEDs used for the 



3.4 Radio to mid-IR flux density ratios 

The correlation of the radio and far-IR luminosities of star 
forming galaxies is well established over five orders of mag- 
nitude (lYun et alj l200lh and is one of the tightest results 
in astrophysics. The continuation of this relationshi p in the 
mid-IR and to z ^ 1 i s also fairly well established l|Garret j 
2002 1 ; lAppleton etafl l2004h . Furthermore we show in the 
Appendix that the total-IR to radio luminosity relationship 
does continue to high redshifts for ULIRGs. This relation- 
ship can be used as another clear discriminant of radio loud 
AGN; all objects with a relative "radio-excess" are very 
likely to be AGN. 

Spitzer provides comparably deep IR observations to 
the radio which allows comparison of the IR properties of 
the sub-mjy radio population. We choose to use the 24 /im 
band of the MIPS instrument as it is more sensitive than the 
70 /im band with a 3er detection limit of 98 /iJy compared 
to 3cr ~ 3 mjy at 70 /jm. This 24 /im detection limit should 
be sufficient to detect all the luminous starburst galaxies 
{L 2 4fj,m > 3 x 10 10 L e ) which we could detect at radio wave- 
lengths out to z ~ 1 and most starburst galaxies at * > 1, 
given the nominal 1.4 GHz detection limit of our parent sam- 
ple. 

In Figure [1] (b) we plot the distribution of mid-IR to 
radio flux density ratios, SW/rni/Si^GHz, against redshift 
for sources in our sample. We compare this distribution to 
redshifted starburst and AGN templates tracks in a similar 
way to section 13.31 On the whole a large number of 24 /mi 
detected sources follow the SFG template tracks suggesting 
that star forming galaxies make up a significant fraction 
of the sub-mjy radio population. There are many obvious 
"radio-excess" and a few "IR-excess" sources which are both 
likely to harbour AGN. The radio-excess sources are likely 
to host AGN due to their radio-loud nature. The IR-excess 
sources are likely to host radio-quiet AGN with hot dust 
dominating in the m id-IR, i.e. radio-q uiet obscured AGN 
such as those found bv lLacv et al. I l|2007r ). but with moderate 
amounts of star formation dominating the radio emission. 
These "IR-excess" sources likely e xplain the high the mid-IR 
to radio flux density ratio seen in iBovle et al.l (|2007f ) due to 
the selection effects of that study. The general trend toward 
lower mid-IR ratios at higher redshifts and the starburst 
SEDs presented sugg est that the low mid-IR ratio found by 
iBeswick et al.l (|2008f ) is due to k-correction effects to the 
observed ratio. 

As in the previous section we use a 5 a cut-off below 
the more extreme starburst. We again get a cutoff value 
of 5 a — 0.6 dex from assuming a median 20% uncertainty 
in both the 1.4 GHz and 24 fim flux density. This cut also 
allows for uncertainty in the SED tracks. 
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Figure 1. (a) 1.4 GHz to near-IR flux density ratio plotted against redshift and (b) 24 /an to 1.4 GHz flux density ratio plotted against 
redshift. Filled symbols have spectroscopic redshifts and open symbols have photometric redshifts. The large blue symbols are those 
radio sources confirmed as AGN from their morphology (triangles), spectral index (squares) or both (stars). We show the location of the 
flux density ratios of the local 6dF-NVSS sample with green dots. The tracks of star forming galaxies with luminosities over the range 
of interest are marked: LIRG (red dotted lines), ULIRG (red dot-dashed lines). See text more discussion of the templates used. W e also 
over-plot the tracks of a radio-loud (blue dash-dot-dot-dot lines) and a radio-quiet QSO (blue dashed lines) from lEIvis et al ] l|l994h . The 
cut-off (solid black lines) in both cases is chosen to be 0.6 dex (= 5 a) above or below the more extreme of two the SFG tracks. Sources 
without known redshifts are plotted at z = 0. 



Table 2. Results of the discrimination of sub-mjy radio sources. 
The raw numbers are presented for each galaxy type as well the 
number of with spectroscopic, photometric or unknown redshifts. 
We also present the number of different sources corrected for in- 
completeness in the radio survey and the percentage of the total 
number of sources. 



Classification 


N 


N sp cc/N phot /N ? 


Ncorr 


Percentage 


AGN 


178 


31/88/59 


244 


35.8% 


SFGs 


269 


131/138/0 


436 


64.0% 


Stars 


1 


1/0/0 


1 


0.1% 



3.5 Results of AGN/SFG discrimination 

Using the above four discriminators we can discriminate the 
AGN from the SFGs on a statistical basis. Although most 
of these discriminators are mainly just positive AGN identi- 
fiers we can be reasonably confident that we have removed 
nearly all the AGN, at least statistically. In fact, 58% of the 
sources flagged as AGN are done so by two or more discrim- 
inators. We find that 178/449 sources show clear indications 
that their radio emission is due to AGN activity. One radio 
source is, surprisingly, identified as a star from optical spec- 
troscopy. These results are summarised in Table [2] along 
with the incompleteness corrected number and percentage. 
This incompleteness factor is a correction for the decrease 
in the sky area probed at faint fluxes due to the decrease in 
sensitivity of our radio map away from the pointing centre. 

We look at the break down of radio AGN classification 
in Table [3] We find that the two most effective methods are 
the flux density ratios which each account for ~ 60 — 80% of 
the sources classified as AGN, but only overlap for ~ 44% 
of the total. Hence, together they account for almost all, 



Table 3. Break-down of the number of radio sources classified as 
AGN by each of our four discrimination methods. The columns 
and rows correspond to the four methods of AGN/SFG discrimi- 
nation used in sections 3.1 to 3.4: radio morphology, radio spectral 
index, mid-IR/radio flux density ratio and radio/near-IR flux den- 
sity ratio. The diagonal represents the number of sources classified 
as AGN by just the corresponding discriminator. The lower left of 
the table, below the diagonal, gives the number of sources flagged 
as AGN for a particular combination of discriminators. 



classification 


morphology 


a 


mid-IR near-IR 


morphology 


44 






a 


12 


32 




mid-IR 


27 


16 


102 


near-IR 


30 


15 


79 137 



90% (160/178), of the radio sources classed as AGN. The 
sensitivity of our K— band and 24 /im observations are gen- 
erally deep enough that we are not biased by non-detections 
in these bands as that would put a source clearly out of 
the SFG regime, except possibly for the mid-IR/radio ratio 
where above redshift 1 we get a few 24 /im non-detections in 
the SFG regime for the very faintest radio sources. The mor- 
phology and spectral index methods are not as powerful (see 
table [TJ, being primarily limited by the sensitivity and cov- 
erage of the current observations. Of the sources classed as 
AGN, these discriminators are together responsible for clas- 
sifying 35.9%, far less than that due to the flux ratio cuts, 
but these two other selection criteria tend to agree with the 
flux ratio methods. 

We also investigated the effect of systematically vary- 
ing the two main discriminators, the flux density ratios plot- 
ted against redshift, by ±ler (= 0.12 dex). When both cuts 
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Figure 2. The 1.4 GHz radio flux density plotted against /—band 
magnitude of all our sources separated into AGN (blue squares) 
and SFGs (red triangles). Radio sources undetected in /—band 
above 3it are marked as upper limits at / ~ 25.9. The lines rep- 
resent the median /—band magnitude for the SFGs (dashed line) 
and AGN (dotted line) as a function of flux density. 



were changed so as to classify more radio sources as SFGs 
we found a swing (i.e. change in AGN/SFG distribution) 
of 4.2% in the number of SFGs from the parent sample and 
when changed so as to increase the number of AGN we found 
a swing of 5.6%. These values are on par with the Poisson 
statistics expected from dividing the parent sample into two 
roughly equal populations (~ 6%). It is possible that some 
or many of the SFGs may contain AGN, but that these AGN 
do not contribute strongly to the radio emission, e.g. the 
"IR-excess" sources in Figure QJb) . It is also possible that 
some radio SFG are really radio AGN, especially as some 
of the discrimination methods do not work so well at faint 
flux densities. The possible number of AGN interlopers is 
difficult to quantify other than varying the discrimination 
criteria (see previous paragraph) without superior data. 



3. 5. 1 Radio sources with faint optical or near-IR 
detections 

Determining the nature of faint radio sources with weak or 
no detections at other wavelengths is naturally difficult. We 
have 60 sources from the parent sample where we are not 
able to determine a redshift, spectroscopic or photometric. 
In fact, 31 of these sources have no detections, > 3 a, at 
any other wavelength. The radio sources with unknown red- 
shifts are included in Figure [T] (both panels) at z = 0, either 
as upper/lower limits or detections. Their positions in these 
plots suggest that they would be mainly classified as AGN 
whatever their redshift, especially in Figure QJa), where our 
AGN/SFG cutoff has only a weak dependence on redshift. 
We considered the possibility that they are low mass SFGs 
at z ~ 1, i.e. not detected in our deep AT— band of IRAC 
data. Such an object would have a SFR of ~ lOOA^yr -1 , 
but a stellar mass < 10 8 5 M o . A few sources with character- 
istics close to these values have been seen, e.g. the hosts of 
GRBs, but the SFR and non-detection in our B— band ob- 
servations implies UV extinction of Av ~ 8 over the whole 
galaxy. While we may have missed a couple of extreme star 



forming sources it is unlikely that many of these sources are 
starbursts at moderate redshift. 

We present the radio flux density verses /—band mag- 
nitude distribution in Figure [2] This plot shows that the 
brightest radio sources (log(5'i.4GHz/M>^J/) > 2.5) are mainly 
AGN, but with a population of star forming galaxies com- 
ing in at lower radio flux densities and brighter optical mag- 
nitudes than the AGN. This result also suggests that the 
optically un-identified and faint sources are most likely to 
be AGN (e.g. at high redshift and obscured). At z ^ 1 
sources would appear compact enough that we should de- 
tect intrinsically bright, but low surface brightness galaxies 
at the depths of our data which have a typical seeing of 
1 arcsec. Below z = 2.2 we expect to detect all sources with 
log(I/2.2pm/i©) > 10, in our near-IR or IRAC data, corre- 
sponding to quite low stellar masses (i.e. log(M/Mo) ~ 10 
depending on the exact mass-to-light ratio) . We plot the me- 
dian /—band magnitude of the SFGs and AGN in Figure [2] 
which both decrease toward fainter radio flux densities. We 
believe the parent sample to be reliable as we would expect 
0.03 of sources to be spurious in pure Gaussian noise, al- 
though the noise in radio maps is not completely Gaussian. 
In S04 we carefully examined visually all the sources to re- 
move those that were clearly due to artifacts from the radio 
reduction. We can't rule out the possibility that one or two 
sources at low S/N and with no optical remain spurious, 
but our results would not change significantly if this were 
the case. 

Figure [2] shows that the optical properties of the par- 
ent sample d o change with radio fl ux density contrary to 
the results of ISimpson et alj l|2006l) who do not see a de- 
crease in the median magnitude with decreasing radio flux 
density. The di f ferenc e between the results here and those of 
ISimpson et all (2006) can mainly be explained by the fact we 
have deeper radio data (by 0.5 dex) and the trend we see is 
strongest at the faintest flux densities. Results at the highest 
flux densities are likely affected by small number statistics. 
Hence, we conclude that most, if not all, of the radio sources 
with optically very faint counterparts are high redshift radio 
galaxies (their optical faintness implying they are obscured, 
type 2 AGN). They are likely to be similar to the classical 
high redshift radio galaxies, but less luminous (i.e. in the 
24 sC logtLi^GHz/WHz- 1 ) sC 26 range ISaiina et alj|2007r i. 



4 RADIO SOURCE COUNTS BY TYPE 

The total Euclidean normalised radio source counts from the 
13 H field were presented in S04 where we found the counts 
to be in agreement with the many other published counts 
at similar flux densities. However, a scatter between the dif- 
ferent counts from different surveys was found which was 
particularly strong around 200 — 300 pJy with the counts 
from the HDF-N l|Richardsl feoOOl ) being significantly low. 
This scatter was attributed to sample variation due to large 
scale stru cture. An independen t reanalysis of the HDF-N 
counts bv lBiggs fc Ivisonl l|2006h revised them upwards, but 
a smaller scatter due to sample variance nevertheless re- 
mains. However, we note that the counts from the 13 H field, 
Phoenix Deep Fiel d Survey (|Hopkins et al.l l2003h. HDF-S 
l|Huvnh et al.ll20o"5h and the three fields from lBiggs fc Ivisonl 
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1.4GHz flux density (mJy) 

Figure 3. The Euclidean normalised 1.4 GHz source counts for 
the 13 H deep survey: total counts (filed circles), SFGs (trian- 
gles) and AGN (open squares). Error bars represent uncertainties 
derived from Poisson statistics and sample variance. We overlay 
AGN and SFG models from S04 and M. Jarvis et al. (in prep. 
J08). Note that a large part of the up-turn in counts above the 
extrapolation of earlier AGN models (S04 etc.) from brighter flux 
densities is explained by SFGs, but that the AGN population 
partly contributes (about one quarter) to the up-turn even in the 
faintest bins. 



20% scatter over the 



from modeling of th e counts rtR< 
Hopkins et all Il998l ; iHopkinsT H 



1993; 



2004; 



(2006) are consistent with only a 10 
70 — 500 /iJy flux density range. 

We can now re-calculate the source counts at 1.4 GHz 
separately by type, AGN and SFG, correcting for incom- 
pleteness of the radio survey as in S04. As in the original 
presentation of these counts we only take sources above 5 a. 
Uncertainties include Poisson statistics and sample variance 
(which we discuss in § 15. 2p . 

The counts by type are presented in Table U and Fig- 
ure[3] We use slightly different binning than in S04 to obtain 
the most uniform distribution of numbers for both AGN and 
SFG populations per bin. We find that there is a rapid rise 
in the contribution to the co unts from SFGs as exp ected 

Rowan- Robinson e t alj 1 
_ .2004 ISevmour et all 
Huvnh et al.ll2005l ). The star forming galaxies are the more 
numerous population below 120 pJy although the AGN still 
contribute approximately one quarter of the total counts 
even in the faintest bin. So while SFGs dominate, AGN do 
have a small contribu tion to the up-turn, b ut maybe not as 
much as suggested bv lSimpson et ail (|2006l ). The fraction of 
SFGs in the two lowest flux density bins, 72 ±22%, is consis- 
te nt with the ~ 70% of < 0.1 mJy sources found to be SFGs 
bv lMuxlow et all (|2005l ) in high resolution (0.2 — 0.5arcsec) 
MERLIN+VLA imaging of the HDF-N by. We note that no 
sources were detected between 27 — 40 fiJy (S/N> 8er) in the 
inner 3.4 x 3.8 arcmin in that work. Those authors inferred 
that the ~ 8 sources expected from the extrapolation of the 
source counts were resolved out on scales > 0.5 arcsec and 
hence, in our opinion, are likely to be mostly SFGs, i.e. a 
SFG fraction close to 100%. 

The sharp disconnect in the number of SFGs from 0.2 
to 0.1 mJy is perhaps due to the fact that our AGN selec- 
tion method may be more efficient at higher flux densities, 
but it most likely represents the sharp rise in the luminos- 



Table 4. Euclidean Normalised 1.4 GHz Source Counts: Total 
and by Galaxy Type. The first column shows the range of flux 
densities of each bin. The second column is the middle of the 
bin in log space. The last three columns indicate the Euclidean 
normalised source counts: total and for AGN and SFGs separately 
for the first time with their associated uncertainties. 



Flux Range 


Flux 




Counts 




(MJy) 


(MJy) 




(sr-ijy 1 - 5 ) 








total 


AGN 


SFGs 


37.5 - 55.7 


45.7 


3.48 ± 0.53 


0.98 ± 0.22 


2.50 ± 0.41 


55.7- 89.6 


70.7 


3.94 ± 0.54 


1.11 ± 0.21 


2.83 ± 0.41 


89.6 - 156.1 


118.3 


4.12 ± 0.42 


1.61 ± 0.24 


2.51 ± 0.32 


156.1 - 294.1 


214.2 


2.73 ± 0.46 


1.82 ± 0.37 


0.91 ± 0.27 


294.1 - 600.0 


420.1 


3.46 ± 0.79 


2.53 ± 0.64 


0.99 ± 0.50 



ity of a typical SFG from z = to z = 1 and the small 
number statistics in the highest flux density bins. The ef- 
fect of slightly varying the flux ratio cuts in the two main 
AGN/SFG discrimination methods (as described in § I3.5|l 
does not qualitatively change the results in Figure [3] When 
changing the cut-offs in flux ratio by 0.12 dex in favour of 
increased SFG numbers the swing in the faintest flux density 
bin is 5.2% and when changing them in favour of the AGN 
the swing in this bin 10.4%. These swings are less than the 
quoted uncertainties on the source counts in this bin, ~ 20%. 

We overlay some models for the AGN and SFG con- 
tribution to the radio source counts. We show the SFG 
and AGN model from S04 and the A GN model courtesy 
of M. Jarvis and R. Wilman, based on Ijarvis fc Rawlinga 
(2004). The S04 AGN model is basically an extrapola- 
tion from higher flux d ensitie s of the mod e ls us ed in 
iDunlop fc Peacock! (|l990l ) and iHopkins et all (|l998l ) and 
hence not necessarily valid for these low flux densities al- 
though they were the best available at the time. The Jarvis 
AGN model matches the observed AGN counts well and in- 
cludes a component of radio quiet QSOs based on the evo- 
lution of the X-ray AGN luminosity function. This model 
appears to follow reasonably well the increase in AGN above 
the extrapolation of older models, like that used in S04, from 
higher flux densities. The S04 SFG model is derived from the 
local luminosity function of SFGs (from [Sadler et al"1l2002l ) 
with luminosity evolution of the form (1 + z) 2 ' 5 to z = 2 (and 
constant thereafter), and does not fit the source counts well 
by over-predicting the SFG counts at the bright end and 
under-predicting them at the faint end. The discrepancy at 
the bright end may be due to the small volume probed at 
low-redshift where the brightest SFGs are likely to lie. 



5 PROPERTIES OF THE RADIO STAR 
FORMING GALAXY POPULATION 

5.1 Star Formation Rates and Stellar Masses 



The conversion of radio luminosity to SFR (Bell 2003) relies 
on the well known radio-IR relation (e.g. lYun et al.ll200ll ). 
While we show in the Appendix that the radio-total IR re- 
lation does indeed hold to high redshifts/high luminosities, 
to derive SFRs from the radio/IR correlation in such situ- 
ations requires two subtle corrections. The first is that in 
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Figure 4. SFRs of the radio SFGs from the 13 H field plotted against redshift (a) and stellar mass (b). The different symbols represent 
different stellar masses ranges as indicated in panel (a). In panel (a), the solid line represents the 1.4GHz detection limit and the dot- 
dashed line indicates the 10 times L* of the radio star forming luminosity function taking a luminosity evolution of Q = 2.5 which we 
use in Section 15.31 The black dots are the star forming galaxies from the local 6dF-NVSS sample. 



determining the local correlation, lYun et all (|200ll ) did not 
apply k-corrections, which become significant for the highest 
luminosity galaxies in their sample. The second issue is that 
some of these galaxies have relatively flat radio spectra, so 
correcting those observed at high redshift back to 1.4 GHz 
using a nominal slope of —0.7 does not provide a good es- 
timate. When nominal k-corrections are applied to the Yun 
et al. sample, the best fit ULIRG ratio of I/60^m/ii.4GHz 
is ~ 140. We have used the observations of Condon et al. 
(1991) to determine a more accurate slope. We find a final 
ratio of L^o/im/ Liaghz = 128 for a galaxy at z ~ 2 ob- 
served at 1.4 GHz and corrected as if it had a slope of —0.7. 
This result would imply that the SFRs estimated above us- 
ing the standard local ratio are under-estimated by about 
0.1 dex, a correction we have applied in our analysis . After 
these corrections and our choice of a iKroupal (|200lT) IMF, 
our conversion from radio luminosity to SFR is 0.07 dex less 
than that of Bell (2003), i.e. for a given radio luminosity we 
assume a SFR 0.84 times that predicted by Bell (2003). 

We plot the SFRs of our SFG sample against redshift 
in Figure HI a) using different symbols for objects of dif- 
ferent stellar masses. Stellar mass estimates are calculated 
for the 13 H sample of SFGs by normalising a M82 SED to 
our IRAC 3.6 /im fluxes (A"— band fo r 8 sour c es no t cov- 
ered by the Spitzer data) and using a iKrou pa (2001) IMF 
to determine a rest-frame H— band luminosities and assum- 
ing a M82 mass-to- (i7— band)-light ratio. Similarly, we de- 
rived stellar masses for the 6df-NVSS sample from their ob- 
served K— band magnitudes assuming the same M82 SED 
and mass-to-light-ratio. Clearly using one mass-to-light ra- 
tio for a sample of galaxies with a range of SFRs, amongst 
other properties, is not ideal, but this approach suffices for 
the current investigation. The possible selection effects of 
this choice are discussed in section [5, 3. II 

We find two striking results. Firstly, there are many 
very high SFR, > 3OOM0yr -1 , sources above z — 1 com- 
pared to the local Universe. The 6dF-NVSS Survey cov- 
ers a large area of the sky, 17% equivalent to ~ 7Gpc 3 to 
z — 0.3, hence is not biased against detecting rare, excep- 
tionally high SFR sources. Therefore, the lack of exception- 



ally high SFR sources locally is not a volume selection effect 
and the presence of very high SFR galaxies at high redshift 
is o bserved in surveys at o t her wavelengths (e.g . in the sub- 



Chapman et~atl 120051 ; iMuxlow et ail 12003 ; iPope et al.l 



120061 ). Secondly, all our high SFR galaxies have very high 
stellar masses. We discuss this result more within the frame- 
work of current galaxy evolution models in Section [5.31 

We can compare the SFR against stellar mass for both 
our sample of SFGs and the 6df-NVSS sample in Figure|4jb) 
and we find that a trend of SFR correlating with stellar mass 
is apparent in both samples (i.e. from the local Universe up 
to z < 3). As radio flux limited surveys will always sam- 
ple the highest SFR objects at any given redshift, this re- 
sult is consistent with the idea that there is an upper limit 
to the SFR of a galaxy o f a given ma s s. We see a simi- 
lar trend in the results of iNoeske et al.l (|2007l ) who find a 
"main sequence" relation between SFR and stellar mass at 
a given redshift with increasing SFR observed at higher stel- 
lar masses. The hypothesis that lower mass galaxies can only 
achieve lower maximum SFRs could be explained by the fact 
that more massive galaxies generally have more baryons in 
gas and dust which can act as raw material for the star 
formation. Potentially some of the very high radio luminosi- 
ties could be over-estimated due to uncertainties in both 
the photometric redshift and k-correction (i.e. radio spec- 
tral index). Our radio/near-IR flux ratio cut does affect the 
number of sources detected with high SFR to mass ratio, at 
high masses/SFR/redshift. However, we have demonstrated 
in Section [3.51 that by varying this cut by 1 a = ±0.12 dex 
does not affect the results our AGN/SFG discrimination 
very much. We may have one or two AGN interlopers, but 
we are likely complete in a statistical sense. 



5.2 Star formation history of the Universe across 

< z < 3 

The star forming galaxies are of great interest as their radio 
luminosity provides a relatively unbiased measure of their 
star formation rate independent of the effects of dust. Using 
this tracer of star formation we can probe the global cos- 
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mic s tar formation rate fe.g. iLillv et al.|[l996 ; Mada u~et al.l 
1 19961 ) . We restricted the sample to all sources with z < 3 as 
we only detect two SFGs above z — 3. 

The star formation rate densities are calculated using 
the standard 1/V m ax method which allows for the change 
in detection limit across the relatively broad redshift bins. 
The redshift bins are chosen to be equal in size in log(l + z) 
space and to have at least 50 sources per bin apart from 
the highest redshift bin which has 13 sources. The redshift 
bins are deliberately left relatively wide to mitigate large 
scale structure/sample variance, Poisson statistics and any 
possible systematics in the uncertainties of the photomet- 
ric redshifts. The total star formation rate density in each 
redshift bin requires two corrections: 

• Radio survey incompleteness: The radio survey is not 
complete down to the 30 /ijy detection limit at 1.4 GHz, 
mainly due to the attenuation of the primary beam of the 
VLA away from the pointing center, but also due to other 
instrumental effects such as band-width smearing. Hence the 
radio sources at low signal-noise, < 100 ^iy, need to be 
corrected by a weighting factor representing the sources not 
detected in our survey. The weighting factor for each source 
is a measure of the sky area sampled by the 1.4 GHz survey 
at the flux density of the source. This correction factor is a 
function of flux density and position away from the pointing 
center, i.e our detection limit is ~ 60 /ijy at the edge of the 
field. This correction is described in full in S04 and we apply 
the individual weightings from that work to each source. 

• Star forming radio sources below the nominal l.J^GHz 
VLA detection limit: As shown in panel (a) of Figure 2] our 
survey primarily probes only the high luminosity end of the 
the SFG luminosity function, especially at high redshift. Any 
estimate of the star formation rate density needs to include 
the contribution from sources further down the luminosity 
function and below the detection limit. Clearly the luminos- 
ity func tion must evolve rapidly given the fits to the source 
counts (|Hopkins et al.lll998l ; ISevmour et"alll2004 etc.) and 
the appearance of very high SFR sources at z > 1. 

The evolution of the SFG luminosity function remains 
slightly uncertain, but most results find a luminosity evolu- 
tion of Q = 2— 3, using the form (l+z)®, with negligible den- 
sity e volution 1 Rowan- Robinson et"al]|l993l; ISev mour et al.l 



2004!; lHopkinsll2004l; iHuvnh et alJl200St IMoss et alj |200'3 k 
Haarsma et al.l (|2000f ) use a different functional form of the 
evolution with more free parameters, but find quantitatively 
the same behaviours for lumi nosity and densi t y evol ution as 
the authors above. However, lHaarsma et all (|200Ch do not 
provide uncertainties for their fitted parameters. As there is 
no consensus on the evolution parameter in the literature, 
or the exact form of the evolution, we take a representative 
value of Q = 2.5±0.5, which is consistent with preliminary 
results from fitting the radio LF (D04), and P = for the 
luminosity and density evolution parame ters respectively. 
We th en use the local radio SFG LF of iMauch fe Sadler! 
(2007). We use these evolution parameters to calculate the 
fraction of the luminosity density below our nominal detec- 
tion limit for each redshift bin. We give the corresponding 
multiplicative correction factors and associated uncertain- 
ties in table [S] 

The uncertainties include the following effects: 
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log(l + z) 

Figure 5. The comoving star formation rate density (SFRD) of 
the Universe as a function of redshift from this work (red trian- 
gles) and other surveys. The horizontal error-bars represent the 
width of the bins and vertical error-bars represent the combined 
uncertainties (see text for details). The SFRD from a selection 
of othe r methods (UV , H a , far-IR etc.) taken from a compila- 
tion in iHopkinj £2004), are marked by gray lines indicating bin 
widths and associated uncertainties. Fro m this sample we have 
highlighted (as blue squares) the results of Haarsma et al.l (|2000l ) 
which are also derived from radio observations . The local (z = 
0.05) value of the radio derived SFRD from ijMauch k, SadleJ 
l2007t) is plotted as a red star with uncertainties smaller than the 
symbol size. 



• Poisson uncertainties: We chose the redshift bin size to 
include at least ~ 50 sources per bin apart from the highest 
redshift bin. 

• Uncertainty in the evolution of the Luminosity Func- 
tion: As mentioned earlier in this section, the amount of lu- 
minosity evolution of the LF remains uncertain so we have 
used a representative value of Q = 2.5 ± 0.5. Hence we in- 
clude this uncertainty when derived the uncertainty of the 
SFRD. The effect of including this factor gets stronger at 
higher redshifts as can be seen in Figure [S] and the LF cor- 
rection factor in Table [5] 

• Sample Variance: Sometimes referred to as Cosmic 
Variance, this effect is really quantifying how representa- 
tive our narrow field is of the whole sky. The source counts 
between different fields do vary by 10 — 20% as discussed 
earlier, but accurately measuring the field to field variation 
between different surveys is tricky due to different methods 
used and in some cases different telescopes to correct for 
instrumental effects that dominate at the faint end. This is- 
sue needs to be addressed by the different teams working on 
deep radio surveys although we previously noted in Section|4] 
that several major surveys are now in reasonable agreement. 
Hence, to allow for sample variance we simply add a further 
20% uncertainty. 

We further investigate the effects of changing the two 
main discriminators (the two flux density ratios as described 
in Section l3.5f) . When changing both these cuts in favour 
of the SFGs we get an increase in the densities, but this is 
barely noticeable in the lowest redshift bin and increases the 
highest redshift bin by just 10%. When changing the cuts 
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Table 5. Star Formation Rate Density at different Epochs. The 
first column represent the size of the redshift bin. The second 
column indicates the middle of the bin in log(l + z). The third 
column gives the number of SFGs in each redshift bin. The fourth 
column indicates the multiplicative correction due to the part of 
the LF below the detection limit using the SFG LF with luminos- 
ity evolution of the form (1 + z) 2 ' 5 + -5 . The fifth column presents 
the comoving star formation rate density derived from the SFGs 
and the uncertainty. 



Redshift 


< z > 


# 


LF cor. 


SFRD 


Range 








log(M Q yr- 1 Mpc _3 ) 


0.10 - 0.52 
0.52 - 1.10 
1.10 - 1.90 
1.90 - 3.00 


0.29 
0.79 
1.47 
2.40 
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87 
54 
15 


, oe+3.5% 
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9 fis +21.0% 
2 - bS -20.0% 

4 60 +48 ' 3% 

g 43+81.6% 
-53.2% 


_1 4 8 +0-09 

-i-07i° :i 2 6 

-0.83t°£ 

-o.92i :2 



the other way, in favour of the AGN, the densities decrease, 
but only by a few percent at the lowest redshift and by 23% 
in the highest redshift bin. Clearly the effect of changing 
the selection criteria is strongest in the highest redshift bin, 
but this result is not unexpected since this bin typically 
contains the lowest S/N radio sources and lowest number of 
sources. The flux ratio discriminators are least good at the 
highest redshifts, and the radio morphology /spectral index 
methods cannot be applied because most objects are too 
faint. Furthermore, these changes are generally smaller than 
the combined uncertainties from the effects listed above and 
hence indicate that the AGN/SFG discrimination criteria 
used in Section [3] are not the largest uncertainty in deriving 
the SFRD at the highest redshift. 

The results of our determination of the comoving star 
formation rate density of the Universe as a function of red- 
shift are presented in Table and Figure [S] Thes e results 
agree well with the multi- wavelength sample from iHopkind 
(2004), converted to a Kroupa IMF, showing the rapid rise 
from the local value to a value over an order of magnitude 
larger at z ^ 1, followed by a flattening above z ~ 1.5. 
We highlight the results of lHaarsma et all (|2000t ). as blue 
squares corrected to the same IMF and cosmology, who also 
derived SFR comoving density from a deep VLA observa- 
tion of several f ields. Their results are consistent with ours 
at all redshifts. lHaarsma et al.l (|2000l ) perform similar cor- 
rections for instrumental effects and for sources below the 
nominal detection limit, but have very different AGN/SFG 
discrimination methods and less sophisticated methods of 
redshift determination (i.e. for sources without spectroscopic 
redshifts those authors use estimates from I or HK 1 magni- 
tudes, or have random assignments). Their uncertainties are 
relatively small as they do not include uncertainties in sam- 
pl e variance. Our resul ts are also consistent with the results 
of|l vison et al. (2007b) who use stacking of different galaxy 
types in the radio to derive SFR densities at different red- 
shifts. The largest uncertainty in our work is that from the 
evolution of the luminosity function, necessary to correct for 
the un-sampled part of the LF. 



5.3 Contribution to star formation rate density 
by galaxy stellar mass 

Recent studies of the star formation history of the Uni- 
verse have not only measured the rapid change in global 
SFR density, but also the demographics of the star forming 
galaxies. Many authors have observed the evolution of the 
distribution of star formation moving high mass galaxies at 
high redshift to lower mass galaxies at lower redshifts (e.g . 



Cowie et al. 1996; Pantcr et al 



Panter et al. 2007; Nocske et al 



12004 lJuneau et all 120051 ; 



2007). This effect has been 



termed "downsizing" and is analogous to (and possibly re- 
lated to) the apparent shift of the peak number density of 
AGN to lower r edshifts for lower X-ray luminosities (e.g. 
lUedaet al.ll2003T l. 

We are able to examine "downsizing" too with our 
unique radio-selected sample of star forming galaxies across 
< z < 2. In this section we chose to look just at the 
contribution from the I/i.4GHz > 10 x L* galaxies where 
we are mostly complete up to z = 2, Figure l^a). The 
ii.4GHz > 10 x L* top end of the luminosity function repre- 
sents the top 28% of the total SFR density at any redshift, 
for the best fitting Mauch and Sadler (2007) LF assuming 
pure luminosity evolution. Hence we can divide each SFRD 
bin into the contr ibution by stellar mass in cluding that of the 
local sample from lMauch fc Sa dler ( 2007|). We only show the 
fractional contribution by luminosity/SFR density in Fig- 
ure [6l but note that the distribution by number density is 
very similar. Hence, number density and luminosity density 
are both good tracers of downsizing. 

While the sample shown in Figure [6] represents about 
one quarter of the star formation occurring at a given cosmic 
epoch, the Li.4GHz > 10 x sources are the most active at 
any redshift. We find that there is a dramatic change in the 
stellar mass of the galaxies contributing to this part of the 
luminosity function. The highest redshift bin is dominated 
by massive galaxies, log(M/M0) > 11.25, although at this 
redshift we are not quite complete to Li.4GHz = 10 x L* 
and hence we mark points in this bin as limits. In the lowest 
redshift bin (from the 6df sample) we can see the dramatic, 
rapid rise in the number of the least massive galaxies. We 
are even able to see the rise and fall of the contribution from 
intermediate stellar mass ranges at successive intermediate 
redshift ranges. 

These results are qualitat i vely similar to oth er recent 
results (e.g. I Juneau et al. 2005; P anter et all 120071 ). but are 
difficult to compare quantitatively as they differ in the se- 
lection details. The differences are largely due to the very 
different survey volumes probed and methods of determining 
the star formation history. Juneau et al. (2005) reach lower 
SFRs determined from optical emission lines from sources 
in Gemini Deep Deep Survey (GDDS) over an area 20 times 
smaller than the work presented here, hence their sample 
includes far fewer of the rarer, extremely high SFR objects 
that we find. lPanter et al l (2007) use the data from the Sloan 
Digital Sky Survey Third Data Release (SDSS DR3 covering 
a far greater area than our survey, ~ 5200 deg 2 ) and model 
the galaxy spectra to determine the individual and global 
star formation histories. 

What Figure |S] truly shows us is the fall of the contri- 
bution by massive galaxies to the global SFR density. We 
cannot rule out that the contribution of the least massive to 
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Figure 6. The contribution by luminosity/SFR density to the 
£i.4GHz > 10 X L* range of the SFG LF. SFGs with different 
mass ranges are indicated on the plot following our earlier con- 
vention. This high end of the luminosity function (in the form 
from Saunders et al. 1990) represents about 28% of the SFRD 
at any particular redshift. The points at log(l + z) = 0.02 are 
from the 6dF-NVSS SFG sample and the other points at higher 
redshift are from the 13 H SFG sample. 



the total SFRD remains approximately constant over this 
redshift range. As the most massive galaxies are capable of 
the highest SFRs they dominate the upper end of the lumi- 
nosity function at the high redshift which we explore. It is 
likely that as these massive galaxies are capable of such high 
SFRs they use up their material to form stars very quickly 
and hence quash their star formation rate. 



5.3.1 Selection effects 

We are selecting our star forming galaxies on a SFR indica- 
tor that is completely unbiased with respect to obscuration 
and with little selection on redshift beyond that implied by 
our radio flux density limit. At any redshift we are natu- 
rally detect just the most luminous, and hence highest SFR 
galaxies. One thing to consider is whether the apparent ob- 
served downsizing is simply an effect of missing low-mass 
galaxies in our high redshift bins. SFGs in our lowest mass 
bin only become un-detectable at z ^ 3 in our K— band or 
IRAC data. Low mass SFGs at these redshifts with SFRs 
high enough to be included in our sample would extremely 
high Specific SFRs (> lO^yr -1 ) and we argue in §3.6.1 
that we have are unlikely to have missed any such extreme 
sources. Hence, the low detection rate of low mass galaxies 
in the high and intermediate redshift bins is not due to our 
detection limits. 

One might argue that the trend of fewer high stellar 
mass galaxies at lower redshift may be a volume selection 
effect. However, we can compare the number of SFGs se- 
lected from their high SFR to that expected purely on mass 
selection. In the lowest redshift bin for the 13 H sample, 
0.1 ^ z ^ 0.4, we would expect to detect ~ 0.2 galaxies 
in the highest mass bin in this volume simply using the 
fitted Schecht er function to the lo cal stellar mass function 
of SFGs from IPanter et al.l |2007) (assuming no evolution 



from z = 0. iFontana et al.ll2006l ). This prediction is consis- 
tent with our result of no galaxies found in this bin. In the 
highest redshift bin we would expect to detect ~ 40 high 
mass galaxies using the same method as before, i.e. assum- 
ing no evolution of t he mass function. However, given that 
I Font ana et al.l l|2006l ) suggest that this part of the mass func- 
tion is not in place until z ^ 1 this estimate of ~ 40 may be 
an upper-limit and hence is consistent with our result where 
we detect 31 high mass galaxies in this redshift bin. These 
rough estimates have to be treated as such as the whole of 
the highest mass bin is above the knee of the stellar mass 
function where the space density changes most rapidly and 
is hardest to measure accurately. However, it is encouraging 
to see that number of massive galaxies varies approximately 
as expected. The ratio of massive to low mass galaxies must, 
of course, decrease toward higher redshifts if there is hierar- 
chical evolution in the mass function, but our selection on 
SFR shows the opposite trend with a higher ratio of high 
mass to low mass galaxies at higher redshift. 

A further factor to consider is the choice of a single 
mass-to-light ratio for the entire SFG sample. Mass-to-light 
ratios can vary by over an order of magnitude from active to 
quiescent galaxies, but most of our sample are clearly quite 
active by any standard, hence any variation in the ratio will 
be much less than this amount. Photometric stellar masses 
used here are also quite model dependent unlike direct dy- 
namical measurements. There is quite a debate concerning 
phot ometric stellar masses in the literature at the moment 
(e.g. lMarastonll2005f ). We note, though, that relative photo- 
metric stellar masses are more accurate than absolute pho- 
tometric stellar masses so the general trends examined in 
this section and the next are genuine cosmological effects. 

5.4 Characteristic Times 

The specific star formation rate (SSFR) of a galaxy is de- 
fined as its star formation rate per unit stellar mass, and 
is representative of how active a galaxy is compared to its 
past properties. For all but the most massive galaxies our 
survey is biased against low SSFRs as we are selecting on 
SFR only. The inverse of the SSF R has been referred to as 
the characteristic time scale (e.g. lJuneau et al. 2005). This 
number has units of time and gives an idea of the nature 
of the star formation going on in the galaxy. For example if 
at a given redshift the characteristic time of a galaxy is less 
than the age of the Universe (or the time since the galaxy 
formed) then the current SFR is higher than the historic 
mean SFR and the galaxy can be said to be in a "burst 
mode" . Conversely if the characteristic time is greater than 
the age of the Universe at a given time then the observed 
SFR is lower than the mean historic SFR and the galaxy 
can be said to be in a "quiescent" mode. In Figure[7]we plot 
the characteristic time as a function of redshift and along 
with the time since the Big Bang, z — 5 and z = 3. Most of 
our star forming galaxies can be classified as being in "burst 
mode", although this classification becomes less certain at 
high redshift as it more strongly depends on the redshift of 
galaxy formation. This result implies that the current ob- 
served star formation rates in most of the galaxies selected 
here are likely due to some trigger event beyond the initial 
gravitational collapse and formation of the galaxy, perhaps 
AGN activity or galaxy-galaxy interaction. The fact that we 
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Figure 7. The "characteristic time" (the inverse of the SSFR) 
of the radio star forming galaxies plotted against redshift. The 
slowly decreasing solid line indicates the age of the Universe at 
a given redshift and the dotted (dot-dashed) line indicates the 
time since z = 5 (3). Galaxies below these lines are forming stars 
at a faster rate than they were in the past and are said to be 
in a "burst mode" and those galaxies above the line are said 
to be in "quiescent mode". For reference, the Milky Way has a 
characteristic time scale of ~ 100 Gyr putting it firmly in the 
quiescent region of the figure (and off the plot), whilst M82 has 
a characteristic time of ~ 0.7Gyrs, well in the local burst region. 



detect the most SFGs in "burst mode" is simply a function 
of the depth of the survey of the parent sample. All but the 
lowest redshift sources have SFR > lOMoyr" 1 , already an 
elevated value for most local galaxies. A much lower radio 
flux density limit is necessary to detect significant numbers 
of "quiescent" galaxies at any kind of distance. 



6 CONCLUSIONS 

We present the analysis of a multi-wavelength follow up to 
one of the deepest radio surveys currently published. These 
data provide spectroscopic and photometric redshifts which 
are vital in discerning the nature of the sub-mjy radio pop- 
ulation. With a philosophy of minimising the assumptions 
about the nature of this population we consider four dis- 
criminators between AGN and SFGs based on purely the 
radio, or relative radio, properties of each source. These 
four discriminators are: radio morphology, radio spectral 
index, radio/near-IR and mid-IR/radio flux density ratios. 
We find that 178 objects in the parent sample were clas- 
sified as AGN by at least one of our AGN discriminators 
and ~ 90% of these are classified by at least one of the 
flux ratio methods methods. Our radio-selected parent sam- 
ple contains a roughly 40/60 mix of AGN and SFGs, but 
with AGN dominating at high radio flux densities and SFGs 
dominating at faint flux densities. By recalculating the Eu- 
clidean normalised source counts by source type we find that 
the up-turn of the counts below 1 mjy is mainly due to 
SFGs, but still with a ~ 25% contribution from AGN at 
the faintest flux densities. This result is consistent with pre- 



vious model fits to the source counts and estimates from 
ultra-deep MERLIN imaging. 

We find many galaxies at high redshift which have SFRs 
exceeding those found in the local Universe and that these 
galaxies tended to be more massive than SFGs at lower red- 
shift. We also find a trend of SFR correlating with stellar 
mass at all redshifts from the 6df-NVSS sample and the 13 H 
sample, i.e. for a galaxy at a given redshift and stellar mass 
there appears to be an upper-limit to the SFR possible. Af- 
ter correcting for radio SFGs below our nominal detection 
limit (4cr = 30 fiJy at 1.4 GHz), we derive the comoving 
star formation rate density as a function of redshift. We find 
our results are consistent with those derived by methods at 
other wavelengths. We are able to look at the population of 
sources in the top end of the luminosity (Li.4ghz > 10 x L*) 
as a function of redshift and find that the typical mass of 
the sources making up the top one quarter of the star form- 
ing luminosity function changes dramatically with redshift. 
We also find that the characteristic time of the SFGs was 
generally low implying an enhanced current SFR for most 
galaxies compared to their past mean SFR. 

The Radio Morphology and Radio Spectral Index dis- 
criminator methods should be just as powerful as the flux ra- 
tio methods if deeper high resolution and multi-wavelength 
data were available. The mid-IR/radio flux density ratio 
would be improved by deeper mid-IR data. Radio luminos- 
ity, not used as a discriminator in this work, will become a 
powerful discriminator when larger volumes are probed in 
deeper and wider surveys. Unless the radio sources unde- 
tected at any other wavelength are extremely obscured (i.e. 
Ay zZ 8), then we are able to reject the hypothesis that they 
are low-mass SFGs at z ~ 1 because none of these radio 
sources have detections in our sensitive it, B, and g' — band 
imaging. Such sources are most likely high redshift (z > 2), 
obscured type-2 AGN. 

These results imply that deeper future radio surveys 
from eVLA, eMERLIN,' LOFAR and (eventually) SKA will 
mainly detect SFGs, not only high SFR objects at extreme 
redshifts, but also low luminosity, more quiescent SFGs at 
the redshifts probed in this work, giving a better view of 
the star formation in the distant Universe from radio data. 
The fraction of AGN is still non-negligible hence discrimina- 
tion between AGN and SFGs will be the principal challenge 
in exploiting future surveys. This discrimination is by no 
means easy and, at least at present, largely statistical. Fu- 
ture surveys will need to be backed up by multi-wavelength 
data, particularly mid-IR and near-IR, photometry, to con- 
fidently separate the SFGs and AGN, as well as deeper low 
frequency and high resolution radio data. 
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Figure Al. (This figure is available as a gif) Observed rest- 
frame SEDs of local star forming ULIRGs (as indicated in figure 
and normalised at 260 /im) compared to far-IR, submm and ra- 
dio observations of submm galaxies (SMGs, marked by dots) at 
1.4 < z ^5 excluding SMGs with likely AGN. This plot indicates 
that for a carefully selected sample of high redshift star forming 
ULIRGs the radio-IR luminosity relation show no obvious devia- 
tion from that found locally. 
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well with the envelope of the SEDs of the local ULIRGs. In 
the radio, there is no apparent offset from the local galax- 
ies, but the scatter is considerably larger than for the local 
ULIRGs. However, a significant part of this larger scatter is 
likely due to measurement errors, since many of the 1.4 GHz 
flux densities have relatively low signal to noise. A further 
explanation for the scatter is that the high redshift star- 
bursts have a wider range of radio spectral indices. 

This paper has been typeset from a TFjX/ I^TpX file prepared 
by the author. 



APPENDIX A: REDSHIFT DEPENDENCE OF 
THE RADIO-IR RELATION 

The radio-IR relation refers to the very close correlation, 
over five orders of magnitude, of t he radio and I R lu- 
minosities of star forming galaxies (jCondon et al.l Il99ll : 
lYun fc Carilil 2002) as the radio and IR luminosities are di- 
rectly related to star formation. To use radio data as a mea- 
sure of star formation, we need to determine if the radio-IR 
relation, from which radio/SFR calibrations are ultimately 
derived, holds or changes above z = 1 (it has been show n 
to hold up to z ~ 1 in the mid-IR. lAppleton et al.ll2004T ). 
In fact, t wo recent studies have indicated that a change is 
possible (|Kovacs et alj|2006l ; Ivlahakis et al.ll2007T I. We pro- 
ceeded somewhat differently from these studies to avoid two 
possible sources of bias. Firstly, we take a sample of radio- 
selected, high redshift galaxies (which have ULIRG or near- 
ULIRG luminosity) and a sample of local ULIRGs only, since 
the relation has not shown to be luminosity- independent 
into the ULIRG range (i.e. local samples from Yun et al. 
2001 and Bell 2003 are restricted to sub- ULIRG sources, 
log(Li.4G.Hz/WHz -1 ) ^ 24). Secondly, we compared the ra- 
dio and IR rest-frame wavelength fluxes between the local 
ULIRGs and high redshift galaxies, since free-free absorp- 
tion can flatten the radio spectrum and make extrapolations 
back to rest-frame 1.4 GHz uncertain. 

First, we assembled a local comparison sample of 
ULIRGs do minated by s tar fo r mation (according t o the 
studies of iFarrah et all 120031 ; lArmus et all 120071 ) and 
with radio measure ments at least at 1.4 and 8.4 GHz 
|Condon et al.l Il99ll ) . There are seven galaxies that meet 
these requirements; Arp 220, IRAS 1056, IRAS 1211, IRAS 
1434, IRAS 1525, IRAS 17208 and IRAS 2249. Figure \M 
shows their spectral energy distributions, normalized at 
260 jiv ci. We over-plotted the measurements of lKovacs et al.l 
(2006) of high redshift galaxies at 1.4 GHz and 350 and 
850/im. We have not shown the measurements of galaxies 
at z < 1.4, and we also rejected measurements of three sys- 
tems that are likely to be influenced by AGN and one more 
not detected at 350 /im: numbers 2, 3, 9, and 10 in their 
Table 1. Where they were available, we used the improved 
1.4 GHz radio flux densities measurements from Biggs & Ivi- 
son (2006) instead of those from Kovacs et al. (2006) . 

We normalised the plotted points to provide a good fit 
at both 350 and 850 fim (observed) to the SEDs of local 
ULIRGs. As Figure 1 shows, the radio and IR flux densities 
indicated for the high redshift galaxies (1.4 ^ z ^ 3.4) agree 
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